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Abstract. This article reports on recent searches for single-top-quark production by the CDF 
collaboration at the Tevatron using a data set that corresponds to an integrated luminosity of 
955 pb~^. Three different analyses techniques are employed, one using likelihood discriminants, 
one neural networks and one matrix elements. The sensitivity to single-top production at the rate 
predicted by the standard model ranges from 2.1 to 2.6 cr. While the first two analyses observe a 
deficit of single-top like events compared to the expectation, the matrix element method observes 
an excess corresponding to a background fluctuation of 2.3 a. The null results of the likelihood and 
neural network analyses translate in upper limits on the cross section of 2.6 pb for the t-channel 
production mode and 3.7 pb for the s-channel mode at the 95% C.L. The matrix element result 
corresponds to a measurement of 2.7lj'3 ^'-'^ combined t- and s-channel single-top cross 
section. 

In addition, CDF has searched for non-standard model production of single-top-quarks via the 
s-channel exchange of a heavy W' boson. No signal of this process is found resulting in lower mass 
limits of 760 GeV/c^ in case the mass of the right-handed neutrino is smaller than the mass of the 
right-handed W' or 790GeV/c^ in the opposite case. 

PACS. 14. 65. Ha Top quarks - 12.15.Ji Applications of electroweak models to specific processes - 
14.70.Fm W bosons - 12.60.Cn Extensions of electroweak gauge sector 



1 Introduction 

In pp collisions at the Tevatron top quarks are mainly 
produced in pairs via the strong force. However, the 
standard model also predicts the production of single 
top-quarks by the weak interaction via the s- or t- 
channel exchange of an off-shell W boson. While early 
Run II searches by the CDF and D0 collaborations, 
based on data sets corresponding to 162 pb^^, 230 pb^^ 
or 695 pb^^ of integrated luminosity, did not find evi- 
dence for single-top production [UlllIS] , one of the lat- 
est D0 analyses [4 using data with Lint = 1 fb"^ ob- 
serves an excess of single-top-like events of 3.4 cr. 

The single-top production cross section is predicted 
to as+t = 2.9 ± 0.4 pb for a top quark mass of 175 
GeV/c^ which is about 40% of the top-antitop pair 
production cross section. The precise measurement of 
the production cross section allows the direct extrac- 
tion of the Cabbibo-Kobayashi-Maskawa matrix ele- 
ment \Vtb\- Moreover, the search for single top also 
probes exotic models beyond the Standard Model. New 
physics, like flavor-changing neutral currents or heavy 
W bosons, could alter the observed production rate [6] . 

The main obstacle in finding single-top is however 
not the production rate of the signal but the large 
background rate. After all selection requirements are 
imposed, the signal to background ratio is approxi- 



mately 1/16. This challenging, background-dominated 
dataset is the main motivation for using multivariate 
techniques. 

Furthermore, single-top analyses are an important 
stepping stone for the Higgs boson search in the WH 
channel, since the two processes have the same experi- 
mental signature. Understanding the background rates 
and background composition in the W-|-jets sample us- 
ing single-top techniques will be a prerequisite for a 
successful WH analysis. 



2 Standard Model Searches 

In this article we present three new CDF searches for 
standard model single-top production and one new 
search for a W boson decaying into tb. All four anal- 
yses are based on the same event selection and use 
the same Run II data set corresponding to an inte- 
grated luminosity of 955 pb~^. The event selection ex- 
ploits the kinematic features of the signal final state, 
which contains a top quark, a bottom quark, and pos- 
sibly additional light quark jets. To reduce multi-jet 
backgrounds, the W originating from the top quark 
is required to have decayed leptonically. One there- 
fore demands a single, isolated high-energy electron or 
muon {Erie) > 20 GeV, or Prifi) > 20 GeV/c) and 
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Table 1. Expected number of signal and background 
events and total number of events observed in 955 pb~^ 
in the CDF single-top dataset. 



Process 


A'" events 


Process 


A'' events 


W + bb 
W + cc 
Wc 

Mistags 


170.9 ± 50.7 

63.5 ± 19.9 

68.6 ± 19.0 
136.1 ± 19.7 


non-W 
tt 

Diboson 
Z + jets 


26.2 ± 15.9 
58.4 ± 13.5 
13.7 ± 1.9 
11.9 ±4.4 


Total background 


549.3 ± 95.2 


t-channel 


22.4 ±3.6 


s-channel 


15.4 ±2.2 


Total prediction 


587.1 ±96.6 


Observed 






644 



large missing transverse energy from the undetected 
neutrino ^x>25 GeV. Electrons are measured in the 
central and in the forward calorimeter, jryj < 2.0. To 
further suppress events in which no real W boson is 
produced, called non-W background, additional cuts 
are applied. The cuts are based on the assumption that 
these events do not produce by nature but due to 
lost or mismeasured jets. Therefore, one would expect 
small and small values of the angle Z\0 between 
and a jet. Wc further reject dilepton events from 
Z decays by requiring the dilepton mass to be outside 
the range: 76 GeV/c^ < M« < 106 GeV/c^. 

The remaining backgrounds belong to the follow- 
ing categories: Wbb, Wcc, Wc, mistags (light quarks 
misidentified as heavy flavor jets), non-M^ and dibo- 
son WW, WZ, and ZZ. We remove a large fraction of 
the backgrounds by demanding exactly two jets with 
Et > 15 GeV and jryj < 2.8 be present in the event. 
At least one of these two jets should be tagged as a 6 
quark jet by using displaced vertex information from 
the silicon vertex detector (SVX) . The numbers of ex- 
pected and observed events are listed in table [TJ 



ing b jet (highest in Et) is included in this distribu- 
tion. The neural net b tagger gives an additional han- 
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2.1 Likelihood Discriminant Analysis 

One multi-variate analysis uses likelihood discriminants 
to combine several variables to a discriminant to sep- 
arate single-top events from background events. One 
likelihood discriminant is defined for the t-channel, one 
for the s-channel search. Seven or six variables are 
used, respectively. The likelihood functions are con- 
structed by first forming histograms of each variable. 
The histograms are produced separately for signal and 
several background processes. The histograms are nor- 
malized such that the sum of their bin contents equals 
1. For one variable the different processes are combined 
by computing the ratio of signal and the sum of the 
background histograms. These ratios are multiplica- 
tively combined to form the likelihood discriminant. 

One of the variables used in the analysis is the 
output of a neural net b tagger. In figure [l}i the dis- 
tribution of this 6 tag variable is shown for the 644 
data events. In case of double-tagged events the lead- 



Fig. 1. (a) Output distribution of the neural net b tagger 
for 644 candidate events in the W^±2 jets bin. Overlayed are 
the fitted components of beauty- like, charm-like and mistag 
templates The yellow error band indicates the statistical 
uncertainties of the fitted sum. (b) The distributions of 
the t-channel likelihood function for CDF data compared 
to the Monte Carlo distribution normalized to the expected 
contributions. A logarithmic scale is used. 

die to reduce the large background components where 
no real b quarks are contained, mistags and charm- 
backgrounds. Both of them amount to about 50% in 
the W + 2 jets data sample even after imposing the 
requirement that one jet is identified by the secondary 
vertex tagger of CDF [7\. 

The i-channel likelihood function is shown in fig- 
ure [lb. The best sensitivity (expected p-value 2.5%) 
is reached by combining the two likelihood discrim- 
inants in a two-dimensional fit where t-chairnel and 
s-channel are considered as one single-top signal (com- 
bined search). The observed data show no indication of 
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a single-top signal and are compatible with a background- 
only hypothesis (p-value 58.5%). The upper limit on 
the combined single-top cross section is found to be 2.7 
pb at the 95% C.L., while the expected limit is 2.9 pb. 



(a) 



The best fit for the cross sections yields at = 0.2_g 2 pb 
and as = 0.lto;ipb. 



2.2 Neural Network Search 



In the second analysis a neural network is used to com- 
bine 23 kinematic or event shape variables to a pow- 
erful discriminant. The output is given between -1 for 
background-like events and +1 for signal-like events. 
The five most important input variables are reported 
by the neural network package. They are: the output 
of the jet-flavor separator neural network, Mi^t, the di- 
jet mass, 77 of the untagged jet times the charge of the 
isolated lepton and the multiplicity of soft jets in the 
event with 8 GeV j Et i 15 GeV. Figure [5^ shows the 
observed data for the combined search compared to 
the expectation in the signal region defined by a neu- 
ral network output between 0.3 and 1. In the combined 
search where the ratio of i-channel and s-channel cross 
sections is fixed to the standard model value a p-value 
of 54.6% is observed, providing no evidence for single- 
top production. The corresponding upper limit on the 
cross section is 2.6 pb at the 95% C.L.. To separate t- 
and s-channel production two additional networks are 
trained and a simultaneous fit to both discriminants is 
performed. The best fit values are at — 0.2tJ;^pb for 
the t-channel and as ~ O.T^J'^pb. The corresponding 
upper limits are 2.6 pb and 3.7 pb, respectively. The 
observed p-value is 21.9%. 



2.3 Matrix Element Analysis 



Another way to discriminate signal from background 
is to compute leading order matrix elements for signal 
and background processes. The measured four-vectors 
of the jets and the charged lepton are used as exper- 
imental input to compute event-by-event probability 
densities. Constraints on energy and momentum con- 
servation are applied. The jet energy measurements 
are corrected to parton level energies using transfer 
functions. To obtain a relative weight for a certain hy- 
pothesis one integrates over jet energies and the mo- 
menta of the incoming partons using the parton den- 
sity functions as integration kernel. The weights for 
the individual hypotheses are combined event-by-event 
by forming a ratio signal and signal-plus-background 
weights. The resulting discriminant is named event 
probability density (EPD). The measured EPD dis- 
tribution is shown in figure [31 compared to the fitted 
event rates for the different processes. A single-top sig- 
nal corresponding to a 2.3 tr excess is observed. The 
associated single-top cross section is 2.7li['3pb. 
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Fig. 2. Single-top search with neural networks at CDF: (a) 
Data compared to the standard model expectation in the 
signal region (neural network outputs larger than 0.3). (b) 
Likelihood contours of the separate neural network search. 



3 Search for a W boson 



Based on a the same event selection as the standard 
model searches CDF has also searched for a W boson 
in the decay channel, W — > tb that appear in mod- 
els with left-right symmetry, extra dimensions. Little 
Higgs, and topcolor 8 . The signal is modeled using the 
event generator pythia [9J. The invariant mass of the 
charged lepton, the reconstructed neutrino and the two 
leading jets is used as a discriminant. The measured 
data are shown in comparison to the standard model 
prediction in figure [3]d. No evidence for a resonant W' 
boson production is found, yielding limits on the pro- 
duction cross section ranging from 2.3 pb at M{W') = 
300 GeV to 0.4 pb at M{W') = 950 GeV. Utilizing the- 
oretical cross section calculations [8 , lower limits on 
the W mass are set: M{W') > 760 GeV if the mass of 
potential right-handed neutrinos is below M{W') and 
M{W') > 790 GeV otherwise. 
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Fig. 3. (a) Event probability density (EPD) based on 
matrix elements. The inset shows the signal region with 
EPD > 0.7. (b) W' tb search: The invariant mass of the 
charged lepton, the reconstructed neutrino and the two 
leading jets. 
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4 Conclusions 



The collaboration has performed searches for standard 
and non-standard model single-top production using 
a data set corresponding to an integrated luminosity 
of 955 pb~^. No evidence for these processes could be 
established. Two standard model searches based on 
likelihood discriminants or neural networks find no ex- 
cess which can be attributed to single-top production, 
while the matrix element analysis finds an excess of 
2.3 cr compatible with single-top production. The over- 
all consistency of all there analysis is only 1%, since the 
analyses feature a correlation between 60% and 70%. 
At present, there are no hints to other causes than 
statistical fluctuations. Recent measurements utilizing 
the matrix element method and likelihood discrimi- 
nant using 1.5 fb~^ show a 3 cr evidence for single top 
production. 

The W search in the single-top channel establishes 
new upper limits of M[W') > 760 GeV or M{W') > 
790 GeV depending on the mass of the right-handed 
neutrino. 



